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Probing the internal rotation of pre-white dwarf stars with 
asteroseismology: the case of PG 0122+200 
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ABSTRACT 

We put asteroseismological constraints on the internal rotation profile of the GW 
Vir (PG1 159-type) star PG 0122+200. To this end we employ a state-of-the-art as- 
teroseismological model for this star and we assess the expected frequency splittings 
induced by rotation adopting a forward approach in which we compare the theoretical 
frequency separations with the observed ones assuming different types of plausible in- 
ternal rotation profiles. We also employ two asteroseismological inversion methods for 
the inversion of the rotation profile of PG 0122+200. We find evidence for differential 
rotation in this star. We demonstrate that the frequency splittings of the rotational 
multiplets exhibited by PG 0122+200 are compatible with a rotation profile in which 
the central regions are spinning about 2.4 times faster than the stellar surface. 
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1 INTRODUCTION 

The upper left corner of the Hertzsprung-Russell diagram 
is populated by a handful of rapidly variable stars, the 
so-called GW Vir stars. They are very hot and luminous, 
hydrogen-deficient pre-white dwarf stars characterized by 
surface layers rich in helium, carbon and oxygen (Werner & 
Herwig 2006) that exhibit nonradial (/(gravity)-modes with 
periods between 5 and 50 min — see, e.g., Winget & Ke- 
pler (2008) and Althaus et al. (2010). In recent years, accu- 
rate asteroseismology of GW Vir stars has started to yield 
details of the internal structure and evolutionary status of 
these stars. On the observational side, the works of Vau- 
clair et al. (2002) on RX J21 17. 1+3412, Fu et al. (2007) on 
PG 0122+200, and Costa et al. (2008) on PG 1159-035 are 
particularly noteworthy. On the theoretical front, important 
progress in the modeling of the internal structure of PG1159 
stars (Althaus et al. 2005; Miller Bertolami & Althaus 2006) 
has made it possible unprecedented asteroseismological in- 
ferences for GW Vir stars (Corsico et al. 2007ab, 2008, 2009). 

Asteroseismology of GW Vir stars provides information 
about the stellar mass, the chemical stratification, the lu- 
minosity and distance, and several other relevant properties 
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such as stellar rotation rate and the presence and strength 
of magnetic fields. Of particular interest in the present in- 
vestigation is the potential of asteroseismology to place con- 
straints on stellar rotation, an important aspect that has 
been proved to be very difficult of assessing by means of 
traditional techniques — mostly spectroscopy. Specifically, 
rotation removes the intrinsic mode degeneracy of a nonra- 
dial g-mode characterized by an harmonic degree £ and a 
radial order k. As a result, each pulsation frequency is split 
into multiplets of 21 + 1 frequencies specified by different 
values of the azimuthal index m, with m = 0, ±1, . . . , it 
(Unno et al. 1989). Rotational splittings in the power spec- 
trum of a compact pulsator were first discovered in the white 
dwarf R 548 (Robinson et al. 1976). Since then, frequency 
splittings induced by rotation have been detected in a num- 
ber of pulsating white dwarf and pre-white dwarf stars. If 
the rate of rotation is slow compared with the pulsation fre- 
quencies, the frequency separation between each component 
of the multiplet is proportional to the rotation velocity of the 
star. This has enabled to derive the mean rotation period of 
a number of white dwarf and pre-white dwarf stars. Inter- 
estingly enough, this approach provides rotation velocities 
much more precise than those inferred from spectroscopy 
(Koester et al. 1998; Kawaler 2004). 

Going one step forward the simple approach described 



2 A. H. Corsico, L. G. Althaus, S. D. Kawaler, M. M. Miller Bertolami, E. Garcia-Berro, & S. 0. Kepler 



Table 1. Properties of PG 0122+200. 



Table 2. The rotational triplets of PG 0122+200. 
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T eff [kK] 

M, [M ] 
logg [cm/s 2 ] 
log(WL ) 
log(H«/-R ) 

M cnv [Mq] 



80±4( a ) 
0.53±0.1< b > 
7.5 ± 0.5< a > 
1.2 + 0.5W 
-1.67±0.26( c ) 



81.54+5 * 

u.bod_ 014 

7 65+ ' 02 
'•DO-0.07 
,,+0.04 
-0.02 
_, 7 n+0.025 

0.019 + 0.006 



1.14"* 



References: (a) Dreizler & Heber (1998); (b) Miller Bertolami & 
Althaus (2006); (c) Estimated from T ff and g, along with M*. 

above, Kawaler et al. (1999) were the first to explore the 
potential of the inversion methods employed in helioseis- 
mology to infer the internal rotation of PG 1159—035, the 
prototype of GW Vir stars. They found that PG 1159-035 
could be rotating slightly faster at the center than at the 
surface, but such a small contrast could be also compatible 
with rigid rotation within the uncertainties of the observed 
splittings. Kawaler et al. (1999) also compared the patterns 
of (m = 0) period spacings with the rotational splittings of 
PG 1159—035 and empirically found that the rotation rate 
of this star must increase with depth. Recently, Charpinet 
et al. (2009) have employed a forward approach aimed at 
determining the internal rotation profile of PG 1159—035. 
They fit the theoretical splittings corresponding to a model 
inspired in the asteroseismological (non-rotating) model of 
Corsico et al. (2008) to the observed frequency splittings 
of PG 1159—035. The forward approach of Charpinet et al. 
(2009) revealed that this star is rotating as a solid body with 
a rotation period of 33.61 ± 0.59 h. 

In this paper we perform a detailed asteroseismologi- 
cal study aimed at placing constraints on the internal ro- 
tation of PG 0122+200, the coolest known GW Vir star 
(T cff = 80 000 ± 4 000 K and logg = 7.5 ± 0.5; Dreizler & 
Heber 1998), using the best existing evolutionary and seis- 
mic models. We present the observations and describe the 
seismological model of PG 0122+200 in Sect. [2] In Sect. [3] 
we explore the internal rotation of PG 0122+200 employing 
the forward approach. In Sect. U we perform rotation inver- 
sions for PG 0122+200 using the Regularized Least Square 
method (Sect. HTTjl . and the function fitting technique (Sect. 
I4.2|l . We conclude in Sect. [5] by summarizing our findings. 



2 THE DATA AND THE MODEL 

PG 0122+200 currently defines the locus of the low- 
luminosity red edge of the GW Vir instability strip. The 
photometric variations of this star were discovered by 
Bond & Grauer (1987). Some observational properties of 
PG 0122+200 are summarized in the second column of Ta- 
ble [1] We employ the high-quality observational data on 
PG 0122+200 gathered by Fu et al. (2007). They have pre- 
sented multi-site photometric observations of PG 0122+200 
obtained with the Whole Earth Telescope (Nather et al. 
1990) in 2001 and 2002. Using their data, together with those 
obtained in previous observational runs, they succeeded in 
detecting a total of 23 frequencies corresponding to modes 
with 1=1 and have unambiguously derived a mean pe- 
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riod spacing of 22.9 s. The 23 frequencies consist of seven 
rotational triplets (m = —1, 0, +1) and two isolated frequen- 
cies with (probably) m = 0. In Table [2] we show the seven 
triplets extracted from table 5 of Fu et al. (2007). The first 
three columns show the radial order, the harmonic degree, 
and the azimuthal quantum number of the eigenmodes. The 
fourth and fifth columns give the pulsation periods and fre- 
quencies, respectively. The sixth column shows the uncer- 
tainty of the frequencies. Finally, in the seventh column the 
observed frequency splittings are listed, while the last col- 
umn gives their uncertainties. The mean value (averaged on 
the seven triplets) of the frequency separations is 3.74 /iHz. 
Adopting this value, Fu et al. (2007) derived a mean rotation 
period of 37.2 hours. 

The theoretical rotational frequency splittings of 
PG 0122+200 were computed adopting the non-rotating as- 
teroseismological model derived by Corsico et al. (2007b) 
on the basis of the modern set of PG1159 fully evolu- 
tionary model sequences computed by Miller Bertolami 
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& Althaus (2006). These models take into account the 
complete evolution of progenitor stars, through the ther- 
mally pulsing asymptotic giant branch phase and born-again 
episode. Corsico et al. (2007b) constrain the stellar mass 
of PG 0122+200 by comparing the observed period spac- 
ing with the asymptotic period spacing and with the aver- 
age of the computed period spacings. Finally, they employ 
the individual observed periods to find a representative seis- 
mological model for PG 0122+200. This asteroseismological 
model reproduces the m = observed periods (see Table [2]) 
with an average of the period differences (theoretical vs. ob- 
served) of < 0.9 s, and represents a substantial improvement 
over those models adopted in previous works. The charac- 
teristics of the asteroseismological model are compared with 
those obtained spectroscopically in Table [T] As can be seen, 
the asteroseismological model accurately reproduces the ob- 
servational data. 



3 THE FORWARD APPROACH: 
ROTATIONAL SPLITTING FITS 

Within this approach the theoretical frequency splittings ob- 
tained varying the assumed rotation profile are compared 
with the observational ones until a best global match is 
found (Charpinet et al. 2009). The goodness of the match 
between theoretical (Su Um ) and observed (Su Mrn ) rotational 
splittings is described using a quality function defined as 
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where we have replaced the subscripts (k, £, m) with a sin- 
gle integer index i that labels the specific splitting (i — 
1, • • • , N s = 14). Each term of the sum is weighted with the 
inverse square of the standard uncertainty (tJi) of the ob- 
served splittings, which are derived from the uncertainties 
in the frequencies given in Fu et al. (2007) and are shown 
in the last column of Table O This is at variance with the 
preliminar study of Corsico & Althaus (2010), in which the 
fits of the rotational splitting were made without weighting 
the terms of the sum, and so, the impact of the different 
uncertainties of the observational data on the final result 
was neglected. The lower the value of \ 2 > the better the 
match between the theoretical and the observed frequency 
splittings. 

The theoretical rotational splittings are computed using 
the expressions resulting from the perturbative theory to 
first order in (the rotation rate) that assumes that the 
pulsating star rotates with a period (P = l/f2) much longer 
than any of its pulsation periods (Unno et al. 1989). Under 
the assumption of rigid rotation (fl constant), the theoretical 
frequency splittings are given by: 



Svkim = -mQ (1 - Cm) , 



(2) 



with m = 0, ±1, . . . , ±£, and Ckt being coefficients that de- 
pend on the eigenfunctions of the pulsation mode obtained 
in the non-rotating case. Such coefficients are computed as 
(Unno et al. 1989): 



Ckt = 



f*'pr 2 [2^ t +g]dr 



(3) 




Figure 1. A propagation diagram of PG 0122+200 showing the 
logarithm of the squared Brunt- Vaisala (solid line) and Lamb 
(dashed line) frequencies. Also depicted is the location of the 
radial nodes of g- modes (blue plus). The loci of the nodes cor- 
responding to the modes exhibited by PG 0122+200 are empha- 
sized with black dots and connected with horizontal lines. 



where £ r and £t are the unperturbed radial and tangential 
eigenfunctions, respectively. In the case of gr-modes, when k 
is large then £ r <C £t, in such a way that Chi — > l/£(£-\- 1) 
(Brickhill 1975). 

If the condition of rigid body rotation is relaxed and 
(spherically symmetric) differential rotation is assumed, £1 = 
f2(r), the frequency splittings are given by (Unno et al. 
1989): 



x T 



Q(r)K M (r)dr, 



(4) 



Kki{r) being the first-order rotation kernels computed from 
the rotationally unperturbed eigenfunctions as (Unno et al. 
1989): 



K kl (r) 



P r 2 {e~^rtt~etW+i)-A} 



(5) 



From Eq. @ it is clear that the frequency splitting for a 
given mode is just a weighted average of the rotation rate 
f2(r) throughout the star, being the rotation kernel Kkt(r) 
precisely the weighting function. 

Note that the perturbative theory to first order in Q, pre- 
dicts symmetric separations of the m/0 components within 
each multiplet with respect to the central one (m = 0) (see 
Eqs. ((2| and Jl])). Therefore, in this work we are neglecting 
the departures from symmetric frequency splitting within 
the triplets centered at II ~ 560 s and II ~ 610 s exhibited 
by PG O122+2OC0. 

To estimate the region of the star probed by the ob- 
served gr-modes, we first examined a propagation diagram 
of the asteroseismological model of PG 0122+200. In Fig. \T\ 
we plot the logarithm of the squared Brunt- Vaisala and the 



1 There exists a number of reasons for these departures (see Vau- 
clair et al. 2002); for instance, the presence of a magnetic field. 
The additional frequency shift due to a magnetic field is depen- 
dent on m 2 , and thus, it could produce an asymmetry in the 
frequency shifts of the m = +1 and m = — 1 components relative 
to the m = component. 
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Figure 2. The normalized rotation kernel Kki(r) for the case of 
k = 12 (left) and k = 24 (right). 



Lamb frequencies, along with the location of the nodes corre- 
sponding to (/-modes (the zeros of the radial eigenfunctions) 
marked with (blue) plus symbols. The nodes associated to 
the eigenmodes exhibited by PG 0122+200 are emphasized 
with black dots. Note that these modes have nodes in the 
region 0.1 < r/R* < 1, implying that they have an oscilla- 
tory character in almost the whole star. We also examined 
the rotational kernels computed from our asteroseismolog- 
ical model for PG 0122+200. In Fig. [2] we show the nor- 
malized Kkt(r) for k = 12 and k = 24, corresponding to 
the shortest and the longest pulsation periods observed in 
PG 0122+200. As can be seen, the rotation kernels have the 
largest amplitudes at the outer regions of the model, but 
have also appreciable amplitudes (up to w 0.3) throughout 
the Ml model of PG 0122+200, implying that the observed 
gr-modes are sensitive to the entire rotation profile. This is 
in contrast to the case of DBV or DAV stars, in which ro- 
tational kernels sample only the outer regions of the star — 
see Kawaler et al. (1999) — because of the larger degeneracy 
of the core. 



3.1 Rigid rotation 

First, we assumed that PG 0122+200 rotates as a rigid body, 
that is, the rotation rate f2 is constant throughout the star. 
We variated the value of fi from 1.15 x 10~ 7 Hz (P ~ 100 
days) to 2.77x 10 -4 Hz [P ~ 1 minutes) and for each value of 
f2 we computed the theoretical frequency splittings (8v^ lrn ) 
by means of Eq. ([2]), where the coefficients Chi were as- 
sessed for each mode through Eq. ([3]) and not by means of 
the asymptotic relation, Cm ~ 1/1(1 + 1). The results of this 
optimization procedure are shown in Fig. that shows the 
X 2 function versus the rotation rate. The best-fit solution 
corresponds to a rotation rate of = 6.915 /iHz. It corre- 
sponds to a rotation period of P = 40.17 hours, in good 
agreement with the approximate value of P = 37.2 hours 
quoted by Fu et al. (2007). Note that the rotation period (of 
the order of 10 5 s) is much longer than the longest pulsation 
period exhibited by PG 0122+200 (~ 600 s), thus justify- 
ing the use of the perturbative theory to a first order in the 
calculation of 8v]; em . 

3.2 Differential rotation 

Here, we lift the assumption of solid-body rotation. Because 
the exploratory nature of this study, we try very simplified 
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Figure 3. The quality function \ 2 m terms of the rotation rate 
Q for the case in which we assume that PG 0122+200 rotates as 
a rigid body. Note the presence of a well defined solution corre- 
sponding to f2 = 6.915 fiHz (arrow). 



functional forms for Q(r). Specifically, we adopt a family of 
linear differential rotation profiles defined as: 

Q(r) = (a - fi c )r + n c , (6) 

where Q s and S7 C are the rotation rates at the stellar sur- 
face and center, respectivelj0. This family of linear profiles 
includes rotation rates that decrease and increase linearly 
with r and also "flat" rotation profiles (when Q s — fi c ) 
that represent the case of rigid rotation already examined in 
Sect. 13. 1] We performed our optimization procedures varying 
the parameters Q s and Q c in the range — 20 piHz. We com- 
puted the theoretical frequency splittings by means of Eq. 
Q, where the rotation kernels are computed by using Eq. 
(0 ■ The results are shown in Fig. [3] where we prefer to plot 
1/x 2 instead of \ 2 1° emphasize the location of the values 
(f2 c , f2 s ) providing good agreement between observed and 
theoretical frequency splittings. The region of good solutions 
(that is, the smallest values of \ 2 ) nas an elongated shape. 
As can be seen, there exists a unique, well-localized best- 
fit solution at (f2 c ,fi s ) = (10.62, 4.41) fiHz, marked with a 
black dot in the plot. This solution is substantially different 
from rigid-rotation, which should fall at some point along 
the green dashed line. The existence of the best fit solution 
suggests that the central regions of PG 0122+200 could be 
rotating more than twice faster than the surface. 

We studied the sensitivity of our result to each of the 
observed rotational triplets. To this end, we have performed 
a simple experiment: we removed one of the observed triplets 
from our computations, keeping the six remainder, and 
searched for a solution. We started with the k = 22 complex, 
that naively could be thought to have the strongest influence 
on our results due to its large frequency separations. Indeed, 
with an average frequency spacing of 4.53 fiHz, it is signifi- 
cantly larger than the other spacings (see Table [2}. We again 
find differential rotation with the core spinning more than 



2 We have also tried two-zone rotation profiles like those used 
by Charpinet et al. (2009) for PG 1159-035, but we were unable 
to derive meaningful properties of the rotation of PG 0122+200, 
which is symptomatic of the inadequacy of the two-zone rotation 
profiles to represent the internal rotation of PG 0122+200. 
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Figure 4. Contour map of the inverse of the quality function \ 2 
showing the goodnesses of the fits in the O s — f2 c plane. We assume 
that PG 0122+200 rotates according to a linear profile. The color 
scale is related to the value of 1/x 2 ■ Light regions are associated to 
the highest values of 1/x 2 , i.e., the best matches between observed 
and theoretical frequency splittings. The location of the best-fit 
solution is shown with a black dot. The green dashed line indicates 
the locii of solutions corresponding to rigid rotation. 



twice faster than the surface, being the variations in the val- 
ues of Q a and Q c negligible when compared with the case in 
which all the seven splittings are taken into account. In fact, 
in this case, we have fl c = 10.75 pHz and fi s = 4.40 /iHz. 
THe similarity of the solutions is due to its large obser- 
vational uncertainty, which, according to Eq. {l}, strongly 
attenuates its impact. A similar experiment in which we re- 
move the k — 24 triplet leads to very similar results. Clearly, 
these two triplets, which suffer from the largest uncertain- 
ties in the list, have no appreciable influence in our results. 
We repeated this experiment with the rest of the splittings. 
We found that the most critical triplets are, in order of de- 
creasing importance, k — 18, 17, 14 and 15. The frequency 
separations for these triplets are accurately known. In par- 
ticular, if we discard the k = 18 triplet, the solution becomes 
Q c — 16.45 /iHz and Sl s = 0.95 £tHz, that is, strong differen- 
tial rotation. On the other hand, if the k = 17 triplet is not 
considered, the solution turns out be compatible with rigid 
rotation, with £7 C = 7.70 /iHz and fi s = 6.55 pHz. Hence, 
our results rely mostly on these two triplets, for which the 
frequency spacings are well determined (although see Sect. 
153) . 

Finally, we examined the quality of the match between 
the observed and theoretical frequency splittings. We find 
that the best fit solution reproduces the observed frequency 
splittings with a mean difference of sa 0.25 piHz, which is 
reduced to ~ 0.1 pHz when we do not consider the difference 
corresponding the k = 22 triplet. 



3.3 Uncertainties 

In order to assess the uncertainties in the derived param- 
eters, we repeated our optimization procedure adding arti- 



ficial, normally (Gaussian) distributed uncertainties to the 
set of 14 observed splittings, with a standard deviation of 
fnoiso = 0.08 /iHz, which is comparable with the best data 
available for pulsating white dwarf and pre- white dwarf stars 
observed with the Whole Earth Telescope (Kawaler et al. 
1999). We performed about 400 realizations of this type. 
The resulting solutions are shown in the left panel of Fig. [5] 
with blue dots. Assuming that the distribution of solutions 
(f2 c , f2 a ) is also Gaussian, we can estimate its dispersion, a. 
We find that the solutions that deviate more from the best- 
fit solution (free of uncertainties, black circle) are located 
at more than 2er from the line that defines the solutions of 
rigid-body rotation (green line). In other words, rigid rota- 
tion can be discarded at a level of confidence of more than 
2cr. We find Q c = 10.62 ± 1.8 (Mz and Q s = 4.41 ± 1.1 /xHz. 
We have repeated the above analysis adopting a deviation of 
Onoiso = 0.12 /xHz for the distribution of uncertainties in the 
frequency splittings (right panel in Fig. ■ These artificial 
uncertainties are by far larger than the average of the un- 
certainties in frequencies quoted by Fu et al. (2007) in their 
table 4. Even in this extreme case, we can discard uniform 
rotation at a level of confidence of more than 1.5a. Thus, 
the conclusion of differential rotation for PG 0122+200 re- 
mains unchanged even when we consider exaggerated un- 
certainties in the measured frequency splittings. Finally, we 
investigated the effects of uncertainties in the asteroseismo- 
logical model on our results, and we found that they are not 
relevant. 



3.4 The drift of the oscillation frequencies 

Vauclair et al. (2011) have published the results of a com- 
prehensive monitoring of seven oscillation frequencies of this 
star. They report changes of these oscillation frequencies 
over time, with much larger amplitudes and shorter time 
scales than those expected by cooling, although the data 
resolution is rather low, of the order of 1 piHz. We focus 
on Table 4 of Vauclair et al. (2011), which shows the fre- 
quency and amplitude variations of the 7 largest amplitude 
modes of PG 0122+200 corresponding to the triplets cen- 
tered at 2224 fiHz {k = 17) and 2493 fiUz (k = 15). These 
two triplets are the only ones (out of seven triplets present 
in the star) which have been well documented to exhibit 
changes over time. 

From Table 4 of Vauclair et al. (2011) we have computed 
the average value of the rotational shifts for each triplet as: 



8v(t n ) — 8v n = 



8v (+) + 



where 5v, 



(+) 



= ^(m=+i) — J / (m=o) an d <5f(_) = V{m=o) — 
i/r m= _ 1 \ at a given epoch t n of observation, where n = 
1, • • • , iVcbs, being A D b s = 6 for the triplet centered at 
2224 fiHz and A b s = 8 for the triplet centered at 2493 /iHz. 
In this way, we are discarding the effects of the observed 
asymmetries in the frequency splittings within both triplets 
of PG 0122+200. This is consistent with the fact that, in 
this work, we are using the the perturbative theory to first 
order in f2 for estimate the frequency splittings, that indeed 
does not account for possible departures of uniformity of the 
splittings within a given multiplet (see Sect. [3}. Next, we es- 
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Figure 5. The best-fit solution without uncertainties (black small circle) and the best-fit solution corresponding to 400 realizations of 
the optimization procedure that include artificial noise (blue dots) with different standard deviations (o- no ; ae ). The green line indicates 
the locii of solutions of rigid rotation, where (Q c = f^s)- 



timated an average value of the splittings over the complete 
set of observations, namely 



n=l 



and the fluctuations around this value A n = 5v n — (5v) 
(n = 1, • ■ • , iVobs)- Finally, we computed the mean value of 
these fluctuations around the average value, 



"Ob! 

< A "> = JT~ E i A -i 



and the variance of the fluctuations 



\ 



iVob 



E 

n=l 



(A„ 



We find for the two triplets (A„) = 0.106 - 0.113 /iHz and 
cr A „ = 0.05 - 0.06 /xHz. 

We check our results of non-rigid rotation by consid- 
ering these possible variations of the rotational splittings. 
We have performed new simulations of our optimization 
procedure, in each of them adding Gaussian noise to the 
splittings of the seven triplets exhibited by the star, with a 
standard deviation oWse- To be consistent, we should adopt 
fnoiso ~ o A n ~ 0.055 piHz in order to estimate the uncer- 
tainties of the frequency splittings due to the observed time 
variations, but we prefer to be somewhat conservative and 
adopt (T no is C ~ (A n ) ~ 0.11 fjJIz, thus overestimating to a 
some extent the impact of the frequency drifts on our re- 



sultsB As expected, the results of our simulations are very 
similar to those shown in the right panel of Fig. [S] indicat- 
ing that rigid body rotation can be discarded at a level of 
confidence of more than ~ 1.5cr. 



4 THE INVERSE PROBLEM 
4.1 The RLS inversion fits 

We also investigated the internal rotation rate of 
PG 0122+200 using the Regularized Least Squares (RLS) 
fitting technique (Kawaler et al. 1999) that has been exten- 
sively applied to the case of the Sun — see Christensen- 
Dalsgaard et al. (1990) and references therein. In this 
method, the internal rotation profile f2(r) is obtained by 
inverting the equation (Jeffrey 

rR, 



5v> 



Q.{r)Ki(r)dr 



(7) 



where 6v® is the ith observed frequency separation (i = 
1, • • • ,M, where M is the number of observed multiplets). 
Here, 6v® corresponds to the average separation within 
a given multiplet. Because the problem is intrinsically ill- 
posed, the inversion of the above equation necessarily must 
be regularized (Jeffrey 1988, Kawaler et al. 1999). To do so 



we minimize: 

M 



i=i 



[8v?-6v?] a +\ [ [£a{r)] s 
Jo 



dr. 



(8) 



3 Note that this is actually an estimation of the magnitude of 
time variations of the frequency separations of PG 0122+200, be- 
cause we do not know the changes that could be experiencing the 
frequencies and frequency separations that have not been moni- 
tored by Vauclair et al. (2011). 
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where the second term of the right hand side is, precisely, 
the regularization term. The form of the regularization term 
determines some additional constraints to the solution. For 
instance, if fl(r) cannot have a steep spatial gradient, then 
£ = d/dr. On the other hand, if f2(r) must be smooth, then 



£ = d 2 /dr 2 . 
numerically: 

R» 



The theoretical splittings 8v i are computed 



Q.(f)Ki(r)dr ; 



N 
3=1 



(9) 



where Qj = fl(rj) and Kij = Ki{rj), being the index j 
(j = 1, • ■ • , TV) associated to the radial mesh point in the 
stellar model on which the kernel Ki(r) and the rotation 
rate Q(r) are evaluated, and Wj = rj+i — rj. 

In the least squares method, the derivative of the sum of 
the squared residuals (the function S) with respect to Qj is 
taken and equated to zero. After some algebra, the following 
matrix equation is derived: 



(K K + AH)f2 



K T Y, 



(10) 



where K is a (TV x M) matrix with elements (K)y = 
WjKij/ai, tt is a (TVx 1) vector with components (f2)j = Slj, 
and T is a (M x 1) vector with components (T)i = Svf /a. 
Finally, H is the (TV x TV") regularization matrix, which 
adopts a tridiagonal form (that is, (H)ij = for \i — j\ > 1) 
or a pentadiagonal structure (that is, (H)y = for \i — j\ > 
2) depending on whether £ ee d/dr or £ ee d 2 /dr 2 . Eq. ()10[) 
constitutes a (TV x TV) system of linear equations that must 
be solved for the unknown rotation velocities that mini- 
mize S. 

We have tested the reliability of our RLS scheme by em- 
ploying this technique on "synthetic" (free of uncertainties) 
frequency splittings generated with the asteroseismologi- 
cal model of PG 0122+200 through the forward approach. 
Specifically, we considered rotational splittings correspond- 
ing to consecutive I = 1 g-modes with k = 1, • ■ • , 40. We 
employed the LU decomposition and also the Gauss- Jordan 
methods (Press et al. 1992) for solve the system. Both meth- 
ods give almost identical results. In all of the cases we have 
examined, the inversions are able to recover the input rota- 
tion profile that we used to compute the synthetic splittings, 
provided that an adequate range of values of the parameter 
A is adopted. 

We have applied the RLS method to infer the inter- 
nal rotation profile of PG 0122+200. We have employed the 
M = 7 averaged 1 = 1 splittings. The regularization ma- 
trix corresponds to the smoothing of the second derivative 
of fl(r). In Fig. [6] we show the inverted rotation profiles for 
PG 0122+200 for several values of A. For very small values 
of A, the inverted profiles exhibit strong variations that lack 
physical meaning. However, as the value of A is increased, 
the inverted solution gradually stabilizes. The resulting ro- 
tation profile (corresponding to A > 10~ 2 ) consists of an al- 
most linearly decreasing rotation rate with f2 c ~ 10.75 /iHz 
and Q, s ~ 4.58 /iHz, in excellent agreement with the results 
of the forward approach. The monotonic linear functional 
form characterizing the inverted rotation profiles should not 
be surprising, since we are forcing f2(r) to have a small value 
of its second derivative at the outset. An analysis of the 
uncertainties similar to that performed for the forward ap- 
proach leads to the conclusion that even with the inclusion 
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Figure 6. Inverted rotation profiles for PG 0122+200 corre- 
sponding to different values of the regularization parameter A. 



of uncertainties in the observed splittings, the rotation of 
PG 0122+200 is faster at the central regions than at the 
surface. Specifically, we found £l c = 10.75 ± 2.4 /iHz and 
O a = 4.58 + 1. 7 /iHz. 



4.2 The function fitting method 

We have also made rotational inversions onto a fixed func- 
tional basis — a method called "function fitting", see 
Kawaler et al. (1999) for details. In this inversion technique, 
an explicit assumption about the functional form of fl(r) is 
made. For instance, if Q{r) is assumed to be a polynomial 



in r of degree K — 1, then fi(r) = 



a k r 



Following 



the least squares method, in which we minimize the sum of 
the squared residuals, we obtain the matrix equation: 



(A T A)a : 



A T T, 



where A is a (M x K) matrix with elements 



JV 

(A)ifc = 2J —WjK lj r k j ' 

3=1 



(11) 



(12) 



a is a (K x 1) vector with components (a)t, = au, and Y is 
a (M x 1) vector with components (T)i = 5u® /<Ti. Eq. fTTJ 
is a (K x K) system of linear equations that must be solved 
for the unknown set of values of a*, that minimize S. 

Specifically, we have explored linear (K = 2) functional 
forms for fi(r), defined by the two parameters <zi = Q c and 
ai = (f2 B — O c ). The optimal values we found for these pa- 
rameters are O c = 10.74 ± 2.9 /iHz and O a = 4.57 ±1.8 /iHz, 
in excellent agreement with the RLS fits and also with the 
forward approach. 



5 SUMMARY AND CONCLUSIONS 

Reliable determinations of the rotation rate of white dwarf 
and pre- white dwarf stars are important, because they can 
provide constraints on the theories of angular momentum 
transport from the core to the envelopes of their progeni- 
tors. These theories predict white dwarf rotation rates that 
are much larger than those spectroscopically found (Koester 
et al. 1998; Berger et al. 2005), and even larger than the ro- 
tation rates inferred from pulsating white dwarfs (Kawaler 
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2004) , unless magnetic torques are taken into account (Suijs 
et al. 2008). These magnetic torques provide the necessary 
spin-down of the cores of the white-dwarf progenitors, thus 
helping to understand the slow rotation of white dwarfs. 

In this work, we have explored the internal rotation of 
the pulsating pre-white dwarf star PG 0122+200 on the ba- 
sis of its observed frequency splittings employing a forward 
approach and also two rotation inversion techniques. The 
three methods employed in this work suggest that the inter- 
nal rotation profile of PG 0122+200 is differential, with the 
central regions rotating more than twice faster than stellar 
surface. Rigid body rotation can be discarded for this star 
at a level of confidence of more than 1.5 and 2 a when we 
consider uncertainties in the observed frequency separations 
of 0.12 and 0.08 /iHz, respectively. By averaging our results 
according to the three method employed, we estimate that 
the core to surface rotation ratio is f2 c /f2 s = 2.4 ± 1.3. 

Up to now, only our Sun (Schou et al. 1998) and two 
pulsating (3 Cephei stars — HD 129929 (Dupret et al. 2004) 
and v Eridani (Pamyatnykh et al. 2004) — are known 
to be rotating differentially with depth. Our results for 
PG 0122+200 can be considered as the first asteroseismic 
evidence of differential rotation with depth for an evolved 
star. Finally, the fact that the internal rotation profile of 
the prototypical star PG 1159—035, which is rotating as a 
rigid body, is apparently very different than that obtained 
here for PG 0122+200 could be indicative that PG1159 stars 
are the result of different evolutionary channels and differ- 
ent mass-loss histories, as suggested by recent theoretical 
evidence (Althaus et al. 2009). Our results call for the need 
of similar analysis for other well studied GW Vir stars. 
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